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This work makes the first ever attempt to understand the influence of the black hole background 
space-time in determining the fundamental properties of the embedded relativistic acoustic geome- 
try. To accomplish such task, the role of the spin angular momentum of the astrophysical black hole 
(the Kerr parameter a - a representative feature of the background black hole metric) in estimating 
the value of the acoustic surface gravity (the representative feature of the corresponding analogue 
space time) has been investigated. Since almost all astrophysical black holes are supposed to posses 
some degree of intrinsic rotation, the influence of the Kerr parameter on classical analogue models is 
very important to understand. General relativistic, axially symmetric, non-self gravitating inflow of 
the hydrodynamic fluid onto a rotating astrophysical black hole has been studied from the dynamical 
systems point of view. The location of the acoustic horizon inside such fluid flow has been identifled 
and the associated acoustic surface gravity has been estimated. The dependence of such surface 
gravity has been studied as a function of the Kerr parameter as well as with other dynamical and 
thermodynamic variables governing the fluid flow under strong gravity. It has been demonstrated 
that for retrograde flow, the surface gravity (and hence the associated analogue Hawking tempera- 
ture) correlates with the black hole spin in general, whereas for the prograde flow, the surface gravity 
co-relates with the black hole spin for slow to moderately rotating holes, but anti-correlates with 
the spin for fast to extremely rotating holes. For certain values of the initial boundary conditions, 
more than one acoustic horizons, namely two black hole type and one white hole type, may form, 
where the surface gravity may become formally infinite at the acoustic white hole. The connection 
between the corresponding analogue Hawking temperature with astrophysically relevant observables 
associated with the spectral signature has been discussed. 



I. ON ANALOGUE GRAVITY PHENOMENA AND ACOUSTIC SURFACE GRAVITY 

Despite the remarkable resemblance in between a black hole and a usual thermodynamic system, black holes never 
radiate within the framework of the classical laws of physics. The introduction of the quantum mechanical effects 
radically changes the scenario - black holes radiate due to the Hawking effects [S^, |3l| ■ Observational manifestation 
of the Hawking radiation for the astrophysical black holes is beyond the scope of the present day's experimental 
techniques. In addition, Hawking quanta may posses trans-Plankian frequencies, and physics beyond the Plank 
scale is yet to be realized within the framework of the contemporary knowledge of the theoretical physics. The 
aforementioned issues had been the prime motivation behind looking for an analogue version of the Hawking radiation 
in the laboratory set up and the formulation of the analogue gravity phenomena was thus introduced by establishing 
the profound similarities in between the propagating perturbation within an inhomogeneous dynamical continuum 
and certain kinematical features of the space time as conceived in the general theory of relativity 0, [l^, [s^ [7l|, [tI, |7j] . 
Contemporary research in this direction has gained widespread popularity since it provides a possibility of verifying 
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certain exotic features of the black hole physics (as manifested in the close proximity of the event horizon, including 
the Hawking radiation) directly some experimentally conceivable physical systems in the laboratory set up. Also since 
the smallest length scale encountered in the analogue gravity models are trans-Bohrian rather than trans-Plankian, 
quantum gravity effects usually significant beyond the plank scale may also show up on considerably larger length 
within the framework of the analogue gravity formalism. 

For the common-most example where the aforementioned perturbation is acoustic in nature, and the spacetime 
describing the fluid flow embedding such perturbation is Minkowskian, the acoustic surface gravity k can be computed 
as f?!]: 



(1) 

Evaluated at the aeoustic horizon 



where Cs is the position dependent sound speed (the velocity of the propagation of the embedded perturbation in 

_d_ 



general), u is the bulk velocity of the fluid normal to the acoustic horizon and -S- represents the derivative taken along 



the direction normal to the acoustic horizon. 

In this context it is to be pointed out that for the flat Minkowskian space time describing the fluid flow under 
the influence of the Newtonian gravity as well as for its relativistic generalization [l|, [l^l , the acoustic perturbation 
propagates along the time like curves. Phonons construct the geodesies of acoustic propagation which are null with 
respect to the acoustic metric (the metric governing the propagation of the perturbation), and generate a null surface 
- the acoustic horizon. Such a horizon forms at the transonic point of the fluid flow. Hence the sonic surface is 
identical with the acoustic horizon, and the supersonic flow resembles the acoustic ergo region. It is evident from 
eq. (H)) that the location of the acoustic/analogue 78j horizon and the values of the speed of the propagation of the 
perturbation, the bulk velocity of the embedding continuum, and their space gradients - all evaluated on such a 
horizon - will provide the complete knowledge of the estimate of the analogue surface gravity. A more comprehensive 
expression for the acoustic surface gravity for a generalized analogue system (including the relativistic generalization) 
may be framed as P, [lol ? ] : 



(U - Cs) 



(1 - C,2) dl^ 



(2) 

Evaluated at the acoustie horizon 



where is the Killing field which is null on the corresponding acoustic horizon and ^ = "^^c!^. The algebraic 
expression corresponding to the ^/x^XtJ. can be evaluated once the space time metric describing the fluid flow as well 
as the propagation of the perturbation in a specified geometry is defined. In this context, eq. ([2]) will further be 
discussed in greater detail in the subsequent sections. 



II. ACCRETING ASTROPHYSICAL BLACK HOLES AS ANALOGUE SYSTEMS 



Conventional works on the classical analogue gravity phenomena focus mainly on physical systems not directly 
subjected to the gravitational force. Gravity like effects rather appear to be an emergent phenomena in such config- 
urations. In such systems, only the analogue Hawking effects may be observed and no source for the conventional 
Hawking radiation may therefore be realized. Quite recently, attempts have been made to study the analogue gravity 
phenomena in astrophysical systems where the strong gravity plays a significant role to influence the relevant fea- 
tures of the acoustic geometry for Newtonian, semi Newtonian and complete general relativistic flows, respectively 
[H, [l^ [111, [H, [H, m, Q, m, [5^ [zl] . Special emphasis has been put to study such effects for accreting astrophysical 
black holes. Transonic accretion onto galactic and extra galactic black holes has recently been shown to be a very 
interesting example of the classical analogue systems found naturally in the Universe [ll, [l^ [2l|, 123,] . Such systems 
are unique in the sense that only for such configuration, both the gravitational as well as the acoustic horizons exist 
simultaneously. This makes the accreting black hole candidates the only analogue systems found in the Universe so 
far where the same flow encounters both kind of horizons. Such systems allow to perform a comprehensive study of 
the properties of space time in the close proximity of both the horizons. 

In this context, spherically symmetric as well as the axisymmetric accretion configuration has been studied. The 
non linear equations describing the stationary, inviscid, hydrodynamic accretion onto astrophysical black holes may 
be tailored to form a first order autonomous dynamical systems [1, l60l ? ? - [6^ . The physical transonic solutions 
for such configuration may formally be realized as critical solutions on the radial Mach number versus radial distance 
phase portrait of the flow. For low angular momentum sub-Keplerian axisymmetric accretion of inviscid hydrodynamic 
fluid, multiple critical points (at most three) may appear and the the integral flow solution may contain stationary 
shock ^ [14|, |22|, (22,, 28., 2^ HS, EO, lil, .4 9. - . 51, .76,] . Any real physical global transonic solution joining the event 
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horizon with the source of the accreting material wih perforce have to pass through a saddle type critical point. 
Conventionally, a sonic point is defined at the radial distance where the radial Mach number becomes unity, i.e., u 
becomes exactly equal to Cs- 

Depending on the equation of state used to describe the flow and the flow geometry, a critical point may or may 
not coincide with a sonic point. For axisymmetric accretion flow governed by the adiabatic equation of state, critical 
points are not isomorphic with the sonic points when the flow configuration is assumed to be in vertical equilibrium 
[21I [2^ . A constant height flow or flow in conical equilibrium produces critical points which are identical with the 
sonic points when the adiabatic equation of state is used [f , _52J. For isothermal flow, critical points and sonic points 
are found to be identical irrespective of the flow geometry used [53 |. Hence for adiabatic axisymmetric accretion 
in vertical equilibrium, the critical surface can not be considered as the acoustic horizon, rather the location of the 
acoustic horizon has to be computed by numerically integrating the differential equations describing the accretion 
solution starting from the critical point - detail procedure for which will further be discussed in greater detail in the 
subsequent sections. 

Usually for the low angular momentum accretion studied in this work, a complete multi-transonic solution may 
be referred to the configuration where two transonic solutions through two different saddle type critical/sonic points 
are connected through a discontinuous shock transition. Such a stationary shock usually formed in between the 
outer saddle type sonic point and the middle centre type sonic point. A comprehensive discussion about such flow 
configuration is available in [2lj & [23 |. 

From the analogue gravity point of view, the aforementioned multi-transonicity may be realized through the 
presence of more than one acoustic horizon of different kind. As obvious that the acoustic horizon may be defined by 
the equation u^ — c^ = 0, for stationary axisymmetric flow configuration the acoustic horizon on the equatorial plane is 
a stationary circle of fixed radius. The measure of such radius depends on the initial boundary condition describing the 
accretion flow. Once the corresponding acoustic metric G^i, is appropriately defined over the stationary background 
metric (^^jy, the discriminant of G^^ as defined by 2? = G^^ — GttG^^ may be used as the marker for categorizing 
various sonic states of the flow, and can as well be used to determine whether the acoustic horizon is of black hole 
or of white hole type. Supersonic (subsonic) flow is characterized by the positivity (negativity) of V [l|, [2l| , and the 
change of sign of V occurs at the acoustic horizon. Hence the sonic points are characterizes by I? < — > V > 
transition whereas V > Q — > V < marks the presence of a stationary shock. Following the classification of Barcelo, 
Liberati, Sonego & Visser (2004), it can be shown that for certain values of the initial boundary conditions describing 
the accretion fiow, two acoustic black holes may form at the two sonic points and an acoustic white hole may form at 
the shock location flanked in between two such sonic points. The acoustic surface gravity becomes formally infinite 
at such white holes because both Cg and u changes discontinuously at the shock location. 



III. CONNECTION BETWEEN THE ACOUSTIC SURFACE GRAVITY AND BLACK HOLE SPIN 

In this work, we intend to study the acoustic geometry of an axially symmetric accretion flow in the Kerr metric for 
a conflguration more complex and astrophycially relevant compared to the idealized disc model of uniform thickness 
as studied in [l[ . Our present work explores the low angular momentum relativistic accretion in vertical equilibrium 
(i.e., where the flow thickness is an analytically deflned function of radial distance, as introduced in "25] using the 
dynamical systems approach). The essential motivation will be to investigate how the spin angular momentum of 
the astrophysical black holes influences the estimation of the acoustic surface gravity for multiple acoustic horizons. 
This is to understand how the properties of the curved space time governing the fluid flow (the Kerr parameter a) 
determines the basic features of the relevant acoustic metric (the corresponding acoustic surface gravity k). Such 
information will allow to probe a physical phenomena of profound significance - the role of the background black hole 
metric (as experienced by the embedding fluid flow) in constructing the perturbative curved manifold determining 
the properties of the space time at the close proximity of the acoustic horizon. Hence the nature of the dependence 
of the value of the acoustic surface gravity at the transonic surface of the accreting fluid on the Kerr parameter will 
manifest the interwining between the black hole metric and the acoustic metric. 

Along with the study of the black hole spin dependence of k, we would also like to investigate whether a physical 
quantity as abstract as the analogue surface gravity may be associated with some observable phenomena in connection 
to the black hole astrophysics. As has already been explained, the low angular momentum accretion flow in the Kerr 
metric may generate two acoustic black holes at its two sonic points and an acoustic white hole at the shock location. 
At the shock, the accreting fluid will be heated up and will get considerably compressed. The amount of compression 
will depend on the strength of the shock, where the shock strength Rm — ^7— is conventionally defined as the ratio of 

the pre- to the post-shock radial Mach number. Eventually, jj^ depends on the initial boundary conditions, including 
the black hole spin. The shock strength Rm is a significant marker of the spectral features of the accreting black hole 
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systems [TT|, [l^ IH, [gI] . Stronger is the shock, prominent is the spectral part indicating the radiation from the post 
shock flow. Hence Rm is a good marker for the observational manifestation in the shock formation for sub-Keplarian 
black hole accretion. In this work we intend to study the co-relation between the value of the acoustic surface gravity 
with the shock strength (which will be different for different values of the initial boundary conditions governing the 
flow, including the Kerr parameter) to associate the concept of acoustic surface gravity with the spectral signature of 
the rotating black holes. 

The plan of the paper is as follows: In next section we will provide a brief account of the transonic black hole 
accretion system from the dynamical systems point of view and the details of such accretion in the Kerr metric will be 
provided. Section|V]wiU be devoted to the calculation of the acoustic surface gravity k. In section |VT] the dependence 
of K on a will clarified in detail, and the corresponding connection with the spectral signature will be described by 
providing the dependence of k on various shock related quantities. Finally in section IVTIl we will conclude. 



IV. BACKGROUND FLUID FLOW CONFIGURATION 



The acoustic geometry corresponding to the stationary solution of transonic, non self-gravitating, axisymmetric, 
inviscid hydrodynamic accretion of weakly rotating sub-Keplarian compressible inhomogeneous fluid onto a spinning 
black hole is considered on the background geometry in the Boyer-Lindquist [l2| co-ordinate. The stationarity and the 
axial symmetry correspond to the two generators of the temporal and the axial isometrics, respectively - leading to the 
fact that the total specific energy and angular momentum remains conserved along the streamline. The particular flow 
model for which the stationary axisymmetric solution of the energy momentum and the baryon number conservation 
equations is sought on the equatorial plane, is described in [22j in great detail. In this section such flow configuration 
will first be be introduced in brief for completeness, then the transonic structure for such configuration will be discussed 
at length. For further clarification regarding the transonic properties of general relativistic axially symmetric black 
hole accretion flow in general, the author may refer to, e.g., P, [25l[5^ ]. 



A. Metric Elements and Related Quantities 

In general any radial distance on the equatorial plane has been scaled in the units of GMbh/c^ and any velocity 
involved has been scaled by the velocity of light in vacuum, c. Mbh is the mass of the black hole considered, and 
G = c = Mbh = 1 is used. — I- signature is used along with a azimuthally Lorentz boosted orthonormal tetrad 
basis co-rotating with the accreting fluid. The energy momentum tensor T^'' = {e + p) v^v'^ + pg^'^ of a perfect fluid 
may be defined where e,p and t;^ are the total mass energy density, isotropic pressure in the rest frame, and the four 
velocity of the accreting fluid, respectively. If p is the corresponding rest mass density, the stationary solution of the 
energy momentum conservation equation T.^"^ — and the baryon number conservation equation (pw^).^; — will 
provide two first integrals of motion, e.g., the conserved specific energy £ (the relativistic Bernoulli's constant) and 
the mass accretion rate M, respectively. The actual expression for £ and M will depend on the specific form of the 
metric (and related co ordinate used), the geometric flow configuration, and the equation of state used to describe 
the accretion fiow. 

As for the space time metric, Boyer-Lindquist line element on the equatorial plane can be expressed as [5^ 

ds'^ = gf.^dx^'dx" = - "^—^dt^ + ^{d(j}- ujdtf + ^dr^ + dz^ (3) 

J\ T ZA 

where 

A = - 2r -I- a^, A = r'' + r^a^ + 2ra'^,uj = 2ar/A (4) 

a being the Kerr parameter related to the black holes spin angular momentum. The required metric elements are: 

/ Alxj'^ r^AX A Aoj 

9rr = -J, gtt = y—^ -J- j ^ 9<t"t> = 9t4> = g<t>t = --^ (5) 

The normal derivative d/drj=rj^d^ — ^ — d/dr. The specific angular momentum A (angular momentum per unit 
mass) and the angular velocity D, can thus be expressed as 

A = -^; n^4 = -^^^^±^, (6) 
vt g4,4, + ^gtcj> 
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We also define 



B = 



(7) 



which will be used in the subsequent section to calculate the value of the acoustic surface gravity. 

The flow is assumed to possess finite radial three velocity u (will be designated as the 'advective velocity') 
on the equatorial plane. Considering v to be the magnitude of the three velocity, u is the component of three 
velocity perpendicular to the set of timelike hypersurfaces {S^} defined by — constant. The normalization criteria 



= — 1 leads to the expression for the temporal component vt of the four velocity 



Ar'^A 



(1 



{A2 - AXarA + A^rZ (4a2 - r2A)} 



1/2 



(8) 



Whereas the advective dynamics is assumed to be confined on the equatorial plane only, the thermodynamic profile 
of the flow is essentially obtained by vertically averaging the thermodynamic quantities over the radius dependent 
vertical thickness 



Hir) 



7+1 



of the flow. 

The adiabatic equation of state of the form 



has been used to describe the flow, where 7 = 
constant volume, respectively. 



(7-l)c? 



{7-(l + c2)}{A2i;2_a2(„^_l)} 



p = K 



(9) 



(10) 



a/cv, the ratio of the specific heat at the constant pressure and at 



B. Conservation Equations and the Critical Point Conditions 



The two first integrals of motion, the conserved specific energy £ of the flow and the mass accretion rate M can be 
computed as 



£ = 



(7-1) 



7-(l + c2) 



1 



1 - 



Ar^A 



A^ - 4AarA + \^r'^{A,a?- - r^A) 

1 u 
M = AttA^ H{r)p—= , 
V 1 — 



(11) 
(12) 



by integrating the stationary part of the energy momentum conservation equation and the continuity equation, 
respectively. The corresponding entropy accretion rate S that is conserved for the shock free polytropic accretion and 
increases discontinuously at the shock (if forms), can be expressed as 



1 \ 

7 



AttA^c 



(7-1) 

7-(l + c2) 



Hir) 



(13) 



The conservation equations for £, M and S may simultaneously be solved to obtain the complete accretion profile 
from the dynamical systems point of view (see, e.g., [2^ for further detail). 

The relationship between the space gradient of the acoustic velocity and that of the dynamical velocity can now 
be computed by differentiating eq. 



dr 



1+7 



4 



2 1 /2 + MV'a\ du 
r 2u \ 1 — v? J dr 



(14) 



whereas du/dr may be obtained by differentiating eq. (fTTj) along with taking help of eq. ([T^ 



du 
dr 



2cl 


V - 1 2 VtOX 


X 


(7 + 1) 


_ A r All) _ 


2 


u 




1 


U^VtCF 


(1-7/2) (7+ 1) (1 - u2)ij 


2ij 



(15) 
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where 



i! = x 



2„,2 



1 dA 



A 



a ^ 2\^vt - 



+ A 



dr 



dr 



(1 - fix) dr 



(16) 



Eq. (jlSp as well as eq. (jl4p can readily be identified with a set of non-linear first order differential equations 
representing autonomous dynamical systems, and their integral solutions will provide phase trajectories on the radial 
Mach number M (where M = u/cg) vs r plane. The critical point condition for these integral solutions may be 
obtained by simultaneously making the numerator and the denominator of eq. (|15p vanish, and the aforementioned 
critical point condition may thus be expressed as 



(7 + 1)^ ' 
2'0 — u^vtc 



1/2 

l(.=.c)' 



l(r=re 



xAr 



2r (r - 1) +4A 



1/2 



(17) 



It is to be noted that eq. ()17|) provides the critical point condition but not the location of the critical point(s). It 
is necessary to solve eq. ()11|) under the critical point condition for a set of initial boundary conditions as defined by 
[£, A, 7, a]. The value of and u, as obtained from eq. (|17p. may be substituted at eq. (|lip to obtain a complicated 
non-polynomial algebraic expression for r — rc, rc being the location of the critical point. A particular set of values of 
[£, A, 7, a] will then provide the numerical solution for expression to obtain the exact value of rc- It is thus important 
to know the astrophycially relevant domain of numerical values corresponding to £,X,j and a. 



C. Astrophysically Relevant Domain of the Initial Boundary Conditions 

£ is scaled by the rest mass energy and includes the rest mass energy itself, hence £ — 1 corresponds to a flow with 
zero thermal energy at infinity, which is obviously not a realistic initial boundary condition to generate the acoustic 
perturbation. Similarly, f < 1 is also not quite a good choice since such configuration with the negative energy 
accretion state requires a mechanism for dissipative extraction of energy to obtain a positive energy solution[80|. 
Any such dissipative mechanism is not preferred to study the acoustic geometry since dissipative terms in the energy 
momentum conservation equation may violate the Lorentzian invariance On the other hand, almost all £ > 1 solutions 
are theoretically allowed. However, large values of £ represents accretion with unrealistically hot flows in astrophysics. 
In particular, £ > 2 corresponds to with extremely large initial thermal energy which is not quite commonly observed 
in accreting black hole candidates. We thus set 1 ^ £ ^ 2. 

A somewhat intuitively obvious range for A for our purpose is < A<2, since A = indicates spherically symmetric 
flow and for A > 2 multi-critical behaviour does not show up in general. 

7=1 corresponds to isothermal accretion where the acoustic perturbation propagates with position independent 
speed. 7 < 1 is not a realistic choice in accretion astrophysics. 7 > 2 corresponds to the superdense matter with 
considerably large magnetic held and a direction dependent anisotropic pressure. The presence of a dynamically 
important magnetic field requires the solution of general relativistic magneto hydrodynamics equations which is 
beyond the scope of the present work. Hence a choice for 1 ;S 7 ^ 2 seems to be appropriate. However, preferred 
bound for realistic black hole accretion is from 7 = 4/3 (ultra-relativistic fiow) to 7 = 5/3 (purely non relativistic 
flow), see, e.g., [26] for further detail. Thus we mainly concentrate on 4/3<7<5/3. 

The domain for a lies clearly in between the values of the Kerr parameters corresponding to the maximally 
rotating black hole for the prograde and the retrograde flow. Hence the obvious choice for a is — l<a<l. The allowed 
domains for the four parameter initial boundary conditions are thus [l ^ £ ^ 2,0 < A<2, 4/3<7<5/3, — l<a<l] . The 
aforementioned four parameters may further be classified intro three different categories, according to the way they 
influence the characteristic properties of the accretion flow. [£, A, 7] characterizes the flow, and not the spacetime 
since the accretion is assumed to be non-self-gravitating. The Kerr parameter a exclusively determines the nature 
of the spacetime and hence can be thought of as some sort of 'inner boundary condition in qualitative sense [sH. 
[£, A] C [f , A,7] determines the dynamical aspects of the flow, whereas 7 determines the thermodynamic properties. 
To follow a holistic approach, one needs to study the variation of the salient features of the acoustic geometry on all 
of these four parameters. 
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D. The Critical Velocity Gradients 

Once the value of Tc is computed for an astrophysically relevant set of [£, A,7,a], the nature of the critical point(s) 
can also be studied to confirm whether it is a saddle type or a centre type critical point (see, e.g., [1^ for an analytical 
scheme developed using the eigenvalue problem to accomplish such classification for axisymmetric accretion in the 
Kerr metric) . The space gradient for the advective flow velocity at the critical point can also be computed by solving 
the following quadratic equation 



|(r=re) ^ I (r = rc 



+ C-0, 



where the respective co-efficients, all evaluated at the critical point Tc, are obtained as 



(18) 



h = 



u{\ — u^) 



S2 = 



2^} 



u{l-v?) 1 - 82 



Tl 



T5 



+ 


) 2S165 


(1 -U2) 


2 ^ 


+ 1' 


^3 


1 


2A2 
a 




, s 






r - 1 


2 

H 
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r 


All; ' 


4-0 


T4 = 


2 


c^Ti - 


'{(7- 


'7 + 1 



P = — TT +'^6, 
7 + 1 

^4 = 62 

(7 - 1 - C3) X 



C = -T-5 

' l-u^_ 
52hxvt 



T2 



2cl 


' 2(1 


-62) 




- a2) ^ _ 





1 

A 



2(r-l)^ 



(Tip 

2 

2 \ X 



Alp dr 

^Idx 
2 dr 



2vtu 



(7+l)(l-«2) 



(7 - 1 - C2) + 



(19) 



Note, however, that all quantities defined in eq. (jl9p can finally be reduced to an algebraic expression in rc with 
real coefficients that are functions of [£, A,7,a]. Hence (du/dr)^.^^ is found to be an algebraic expression in rc with 
constant co efficients those are non linear functions of ,A,7,a]. Once r^ is known for a set of values of [£,A,7,a], 
the critical slope, i.e., the space gradient for u at rc for the advective velocity can be computed as a pure number, 
which may either be a real (for transonic accretion solution to exist) or an imaginary (no transonic solution may be 
found) number. The critical advective velocity gradient for accretion solution may be computed as 



'd^l. 



2a 



± - 4aC 



(20) 



by taking the positive sign. The negative sign corresponds to the outflow/self-wind solution on which we would 
not like to concentrate in this work. The critical acoustic velocity gradient {dcg/dr)^.^^ can also be computed by 
substituting the value of in eq. (ITil) and by evaluating other quantities in eq. at rc- 



E. On Transonicity and Multi-transonic Solutions 



As mentioned before, the critical point can be computed by putting the critical point conditions in the expression on 
£ as expressed in eq. ([TT|). The number of critical points obtained is either one (saddle type), or three (one centre type 
flanked by two saddle type) depending on the particular value of [£, A, 7, ajused. Certain [£, A, 7, a]ijicC[f , A, 7, a] thus 
provides the multi criticality in accretion (as well as in outflow) solutions, where the subscript 'mc' stands for multi 
critical. It is, however, to be noted that the radial Mach number at the critical point is a functional of [£, A,7,a]and 
is less than unity, which is obvious from eq. (jl7p . Following our previous discussions, this is the consequence of the 
choice of the geometric conflguration of the accretion flow and the equation of state to study such flow. As mentioned 
earlier, the acoustic horizon is defined as a time-like hypersurface defined by the equation 



(21) 
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The acoustic horizon are thus the collection of the 'sonic' points where the radial Mach number becomes unity. Since 
critical points are not topologically isomorphic with the sonic points in general, acoustic horizon does not form at 
the critical point, and neither [du/dr]^^^ nor {dcs/dr)^.^^. can be used to evaluate the value of the acoustic surface 
gravity. Had it been the case that the critical points would be identical with the sonic points, the acoustic surface 
gravity could easily be evaluated by taking Vc to be r^, the radius of the acoustic horizon, and by directly taking 
the value of Cs,u, {du/dr)^^^ and (dcs/dr)^^^ as evaluated at rh- One thus understands that rh is actually a sonic 
point located on the combined integral solution of eq. and eq. (fT5)) . For inviscid flow, a physically acceptable 
transonic solution can be realized to pass through a saddle type sonic point, resulting the hypothesis that every saddle 
type critical point is accompanied by its sonic point but no centre type critical point has its sonic counterpart. For 
axisymmetric flow in the Kerr metric, a multi-critical flow is thus a theoretical abstraction where three critical points 
are obtained as a mathematical solution of the energy conservation equation (through the critical point condition), 
whereas a multi-transonic flow is a practically realizable configuration where accretion solution passes through two 
different saddle type sonic points. One should, however, note that a smooth accretion solution can never encounter 
more than one regular sonic points, hence no continuous transonic solution exists which passes through two different 
acoustic horizons. The only way the multi transonicity could be realized as a combination of two different otherwise 
smooth solution passing through two different saddle type critical (and hence sonic) points and are connected to each 
other through a discontinuous shock transition. Such a shock has to be stationary and will be located in between two 
sonic points. For certain [£, A, 7, a]nssC[f , A, 7, ajmc, where 'nss' stands for no shock solution, three critical points (two 
saddle embracing a centre one) are routinely obtained but no stationary shock forms. Hence no multi transonicity 
is observed even if the flow is multi-critical, and real physical accretion solution can have access to only one saddle 
type critical points (the outer one) out of the two. Thus multi critical accretion and multi transonic accretion are not 
topologically isomorphic in general. A true multi-transonic flow can only be realized for [£", A, 7, a]ssC[f , A, 7, a]„ic, 
where 'ss' stands for 'shock solution', if the following criteria for forming a standing shock (the relativistic Rankine- 
Hugoniot condition, see, e.g., [U, for further detail including the formulation) gets satisfied 

[[puTu]] = 

[[T,.rvl]^[[{p + eWTl+p]]=0 (22) 

where r„ = f /Vl — is the Lorentz factor. 

The overall procedure to compute the value of the analogue surface gravity, is, however, a bit involved. For a 
set of [£, A,7,a], rc,u\(^r=ro)^'^s\^^^^^y idu/dr)^^^.^ , {dcs/dr)^.^^.^ can be calculated. These values are then used as the 
initial values to find the integral solution for eq. ([T4l - [T5|) . Along the integral solution, the radial Mach number M 
is computed at every r until one obtains the specific value of r for which M becomes exactly equal to unity. rM=i 
is thus equal to Vh- We then calculate u,Cs,du/ dr and dcg/dr at rh and calculate the value of the acoustic surface 
gravity k using such values. The exact expression of k will be derived in the next section. 



V. ACOUSTIC SURFACE GRAVITY IN TERMS OF [£,A,7,a] 

For non-dissipative, barotropic, relativistic irrotational flow, the wave equation describing the propagation of the 
acoustic perturbation can be described as (loj 



-.da 



\G \G 



d^ip = 



(23) 



by linearizing the relativistic Euler equation and the continuity equation around some steady background, ip in eq. 
P3p represents the low amplitude linear perturbation of the velocity potential around the aforementioned steady 
background. The acoustic metric tensor and its inverse may be defined as 



G 



P 

he. 



[g^, + {1- cl)v^v,] ; G'^" 



hcs 
P 



f 



+ {l~-)v^v 



(24) 



where h is the relativistic enthalpy and 3^,^ is the stationary background metric. 

For the steady axisymmetric fiuid fiow considered in this work, the acoustic metric is assumed to be stationary 
and any displacement along the projection of the flow velocity on E„ is assumed to be an isometry, where Ei, is the 
hypersurface of constant v. The acoustic ergo- region in such configuration may be defined as the region where the 
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stationary Killing vector becomes spacelike, i.e., any supersonic region is an ergo region. The stationary limit surface 
is realized as the boundary of the ergo region as defined by the equation 



-ct=0 



(25) 



In terms of the acoustic metric tensor, is defined by the condition Gu 
defined as a timelike hypersurface satisfying the criteria 



0. The acoustic horizon can now be 



(26) 



where rj'^ is the unit normal to the horizon. 

For axially symmetric accretion flow, the acoustic horizon does not co incide with the stationary limit surface in 
general, except where the three velocity component perpendicular to the horizon is considered to be the flow velocity 
(as measured by a co rotating observer) of interest. One thus identifies the advective velocity u on the equatorial 

plane with ^{rj^Vf^) / \J (77^?;^)^ + rji^rif^ , and the acoustic horizon can essentially be defined by eq. (PT|) as explained 

earlier. This further confirms that the sonic horizon (the collection of all the points where M — I) is the acoustic 
horizon. 

One now constructs a Killing vector x'' = + ^4''^ where the Killing vectors and (p^ are the two generators 
of the temporal (constant £) and axial (axisymmetric flow) isometry groups, respectively. When is computed at 
r/j, becomes null on the transonic surface. The norm of the Killing vector Xn may be computed as 



With this, eq. © takes the form 



>^gt<i> 



(27) 



2r + a? 



2\gt4, + X^gtt) 



(1 - 0^2) (r3 + a?r + 2a? - 2\a) 



du 
dr 



dcs 
dr 



(28) 



A knowledge of u, Cg, du/ dr and dcs/dr as evaluated at the sonic point is thus sufficient to calculate k, for a fixed set 
of values of [£, A, 7, a]. 

We now illustrate the procedure to calculate k, for two typical flow configuration - for a montransonic accretion 
solution passing through the inner type sonic point and a multi-transonic accretion solution with stationary shock. 
In figure 1, ABC is the integral solution for transonic accretion obtained by simultaneously solving eq. (|14l - [T5|) . 
and DEF is the integral solution representing the 'self wind'. The radial Mach number is plotted along the Y axis 
and the logarithmic radial distance measured from the black hole event horizon (along the equatorial plane) has been 
plotted along the X axis. The phase portrait is obtained for [£ = 1.2, A = 2.0, 7 = 1.6, a = 0.4]. The intersection of 
ABC and DEF is essentially the critical point Tc located at 3.3276 in units of GM jj-^Qj^ j (? . The critical point r^ 
is obtained by putting the critical point condition as obtained from eq. (jl7p in the energy conservation condition as 
stated in eq. (jlll) . and by solving it numerically for r^- Such r^ and its corresponding functions are then substituted 
back to eq. p7|) to obtain the critical velocity and the critical sound speed. The critical space gradient for u and Cs 
are then obtained from eq. (j20p . [u, Cs, du/ dr, dcs/drj^.^^ is then used as initial value to obtain the integral solution 
of eq. (ITil - [T5]) . Mach number is calculated at each point along the integral accretion solution ABC and the value 
of r for which Af = 1 is identified as r^ {ru < rc for accretion for obvious reason), [u, Cs, du/dr, dcs/dr]^^^^ is then 
computed, and their values are substituted at eq. p8)) to obtain the surface gravity in terms of [£, A, 7, a]. 

Similar exercises may be performed for the integral transonic solution for wind (branch DEF). However, in this 
work we are interested to study the acoustic geometry for the inward moving accretion fiow only, k can also be 
computed for the monotransonic integral accretion solution passing through the outer type sonic point. The acoustic 
surface gravity is way much higher for solutions passing through the inner type sonic point compared to that passing 
through the outer type sonic point[82|. By varying any of the four parameters £,A,7 and a, one can study the 
dependence of the acoustic surface gravity on the flow properties and on the property of the back ground space time 
itself. 

In flgure 2, a multi-transonic fiow with stationary shock has been depicted for 
[£ = 1.00001, A — 2.6,7 — 1-43, a — 0.4]. ABC is the integral accretion solution passing through the outer critical 
point r°"' and DEGF is the lower half of the homoclinic orbit representing the transonic solution passing through 
the inner critical point r™. The vertical line BG represents the shock transition. The shock location is obtained by 
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r 



FIG. 1: Phase portrait corresponding to the mono-transonic accretion ABC (and its associated self-wind DEF) characterized 
by [E = 1.2, A = 2.0, 7 = 1.6, a — 0.4]. The inner type critical point is formed at r — 3.3276, while the sonic horizon is located 
ar r = 2.9503, both in units of GAIbh/c^ ■ See text for further details. 



solving the relativistic Rankine Hugoniot condition as stated in eq. (|22p . CBGED is the combined multi-transonic 
integral solutions with stationary shock. Starting from [u,Cs,du/dr,dcs/dr\^^^ai-Lt and [u,Cs,du/dr,dcs/dr\^^^in, one 
calculates [u,Cs,du/ dr^dcs/ dr\^^^a^t and [u,Cs,du/dr,dcs/dr\^^^iri, respectively. The corresponding surface gravity 

h h 

Kout and Kin are then computed for the outer and the inner acoustic horizons r^"' and r^f^ , respectively. The ratio of 
the surface gravities at the inner acoustic horizon to that of the outer acoustic horizon is then computed as 



At the shock, the pre and the post shock Mach numbers 7\f_ and are the Mach numbers calculated at the 
shock location for the accretion solutions passing through the inner and the outer sonic points, respectively, and 
the shock strength Rm — M^/M^ can be computed for a fixed set of [£,A,7,a]. By varying any one of the four 
parameters f , A,7 or a, the dependence of the ratio of two acoustic surface gravities on the shock strength can the 
studied by plotting k^o as a function of Rm- Similarly, the dependence of k^o on £, A, 7 and a can be studied for the 
multi-transonic configuration. 



VI. DEPENDENCE OF k ON [f. A, 7, a] AND ON VARIOUS SHOCK RELATED VARIABLES 

In this section we will illustrate how the properties of the background space time (the properties of the black hole 
metric manifested apparently through the Kerr parameter a) influences the properties of the perturbed manifold (the 
properties of the acoustic metric as manifested through the measure of the acoustic surface gravity k, or the ratio Kio 
of such k's at the inner and the outer acoustic horizon, respectively) for transonic accretion with and without shock. 
We will also investigate how the dynamical (manifested through the dependence of k on , A]) and the thermodynamic 
(manifested through the dependence of k on 7)propcrties of the fluid flow influences the characteristic features of the 
acoustic geometry. We will consider both shocked and non shocked accretion of the prograde as well as the retrograde 
flow. 

As mentioned earlier, we would like to investigate the variation of k with each of the parameters [£, X,j,a]. It is 
to be understood that while studying the dependence of k on any one of the aforementioned four parameters, say £, 
A, 7 and a has to be kept fixed for the entire range of £ for which the value of k has to be computed. Since for a limited 
non linear range of [£, A, 7, a]one obtains a multi-transonic flow with shock, a continuous range of all parameters can 
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FIG. 2: Phase portrait for the multi-transonic shocked accretion characterized by [£ — 1.00001, A = 2.6, 7 = 1.43, a = 0.4]. The 
vertical line BG demonstrates the shock transition. The inner, middle and the outer critical points are located ai r — 5.2398, 
r = 18.9471 and r = 10829.8325, respectively, in units of GMbh/c^. The inner and the outer sonic horizon are located at 
r — 4.6775 and r = 8449.4640, respectively. The shock is located at r = 72.01. The middle critical point is indicated by the 
symbol '+'. 

not be explored for such flow. Whereas for the mono transonic flow, a wide range of parameters, sometimes even the 
entire range of the allowed values corresponding to such parameters, may be explored. We will illustrate this issue in 
the subsequent section in greater detail. 

A. Black Hole Spin Dependence of k 

Astrophysical black holes can posses a wide range of spin angular momentum [45l - l47i l66l . l69l . pTZj . We would thus 
like to cover a sufficiently large domain of a, from slowly rotating to the near extremally rotating Kerr black holes, 
while studying the spin dependence of the acoustic surface gravity for the prograde flow. However, multi-transonic 
flow characterized by a single set of [£, A, 7] can not cover the entire stretch of the aforementioned domain since the 
shock does not form (i.e., eq. does not get satisfied) for all values of a for a single set of ,A,7]. In this work 
we explore three different ranges of the Kerr parameters for the multi-transonic prograde flow. Flow with different 
ranges of a as are characterized by three different A (=2.6, 2.17 and 2.01 for the three different ranges of a), for the 
uppermost, middle and the lower panel, respectively) for flxed set of values of [£ = 1.00001,7 = 1.43]. 

From recent theoretical and observational flndings, the relevance of the counter rotating accretion in black hole 
astrophysics is being increasingly evident [23,[53,[23- It is thus instructive to study whether the characteristic features 
of the acoustic geometry remains invariant for a direct 'spin flip'. In other words, whether the k - a profile changes 
significantly when the initial boundary condition is altered from [f, A,7,a]to [£ , X,j, —a]. It is obvious that such 
investigation can not be performed for multi-transonic accretion since shock condition can never be satisfied for a 
certain set of [£,A,7, a]as well as for [£,A,7,— a], i.e., for flows characterized by exactly the same value of [i5,A,7,] 
and magnitude wise same but direction wise different values of a. This is because the region of the four dimensional 
parameter space spanned by [£, A, 7, ajresponsible for the shock formation does not allow any parameter degeneracy. 
Monotransonic accretion solutions, however, is not constrained by this issue, and one can obtain such solutions when 
[£, A, 7, ajgets directly altered to [£, A, 7, —a]. However, dependence of k on a for multi-transonic shocked retrograde 
flow has also been performed for in this work, see subsequent sections for further details. 
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1. Monotransonic Accretion 



4.5 I 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




1 I ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

a 

FIG. 3: The location of the acoustic horizon Vh (dashed hne), the corresponding critical point Vc (solid line) and r+ (dot- 
ted line), all in units of GMbh/c^ , as a function of the black hole spin for mono-transonic accretion characterized by 
[£ = 1.2, A = 2.0,7 = 1-6] • 

We first calculate the acoustic surface gravity for mono-transonic accretion corresponding to the phase portrait as 
shown in figure 1. [£ — 1.2, A 2.0, 7 = 1.6] has been used to calculate k for all values of the Kerr parameter ranges 
from a = — ltoa = lto explore the entire range of the prograde as well as the retrograde flow. In figure 3, the 
location of the acoustic horizon rh is plotted as a function of the a. Such variation is shown by the dashed curve. In 
the same figure the variation of Vc with a is shown by the solid curve. As a reference, location of r+ = 1 + Vl — 
as a function of a has also been plotted in the same figure using the dotted curve. The location of the acoustic 
horizon anti-corrclatcs with the black hole spin. Tc anti-correlates with a as well. We define Arjf = (rc — r^) as 
the difference of the location of the critical point and the sonic point (acoustic horizon), respectively, for a fixed set 

of value of [f. A, 7]. Arjf''^''^' non linearly anti correlates with a. For very large value of the black hole spin for co 
rotating accretion, the critical point almost co incides with the sonic point. However, the critical point and the sonic 
point are never identical {Ar^cs'^''^^^0 always) for any value of a. 

In figure 4, we plot the surface gravity as a function of the black hole spin. For retrograde flow, k co relates with 
the black hole spin. For prograde flow, the surface gravity initially increases non linearly with the black hole spin and 
attains a maximum value for a moderately high value of a, after which it falls off non linearly as a is increased further. 

The location of the peak of the 'k — a' graph on the abscissa, i.e., the value of a for depends on the 

choice of [£,X,j]. We find that Om^x'''^ nonlinear anti-correlates with [f. A, 7] whereas Km^x'''' non-linearly correlates 
with [5, A, 7]. For the mono-transonic flow k in general co-rclatcs with [5, A, 7]. This is because anti-corrclatcs 
with [£, A,7], and closer the acoustic horizon forms to the actual black hole event horizon, higher is the value of the 
associated surface gravity. It is interesting to note that similar non-monotonic behaviour of 'k — a' dependence for the 
prograde flow is observed for multi-transonic shocked accretion as well. The surface gravity computed at the inner 
acoustic horizon for such flow exhibits a maximum for a certain high value of the black hole spin. Such results are 
presented in the subsequent sections. 
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FIG. 4: The surface gravity ft evaluated at the acoustic horizon as a function of black hole spin, o, for mono-transonic accretion 
characterized by a fixed set of {£ = 1.2, A = 2.0, 7 = 1.6]. 



2. Multi-transonic Shocked Accretion 

The typical phase diagram for a representative flow topology has already been shown in figure 2. For prograde 
flow, three different ranges of the black hole spin has been studied for same values of [£ = 1.00001,7 = 1.43] but for 
three diflferent values of A - A = 2.6 for a ranging from [a = 0.21067333 to a = 0.47828004], A = 2.17 for a rang- 
ing from [a = 0.85614997 to a = 0.92420954], and A = 2.01 for a ranging from [a = 0.96821904 to a = 0.98999715], 
respectively. For the retrograde flow, [£ = 1.00001, A = 3.3,7 = 1-4] has been used to study the range of a from 
[a = —0.21 to a = —.661200]. Note that the values of [£, A, 7] used here are just some representative set of values for 
which the shock forms for a substantially large range of a. Any other set of values of [E, A, 7] for which a reason- 
able range of Kerr parameter satisfies the rclativistic Rankine Hugoniot condition can be used as well to study the 
dependence of the acoustic surface gravity on the black hole spin. 

In figure 5, we plot the value of the location of the outer acoustic horizon (uppermost panel), the shock location 
(mid panel) and the inner acoustic horizon (lowermost panel) as a function of the Kerr parameter a for prograde 
flow characterized by [£ = 1.00001, A = 2.6,7 = 1-43]. The location of the inner and the outer acoustic horizon non 
linearly anti-correlates with a, whereas the shock location non linearly co relates with a. The variation of the outer 
acoustic horizon is much insensitive on a, which is probably expected since the outer horizon forms at a very large 
distance and a being the inner boundary condition, does not play any significant role in influencing the variation of 
any accretion related quantity at that length scale. On the other hand, variation of the inner acoustic horizon as well 
as of the shock location are considerably sensitive to the black hole spin. Similar variation of the location of the inner 
acoustic horizon r^", the outer acoustic horizon r^"* and the shock location Vsh on the black hole spin can be studied 
for other set of ranges of a for both prograde as well as for the retrograde flow. The overall ' [rj^", Ts/i, r^"*] — a' profile 
remains more or less the same for any range of a used - although the numerical values of [r™, rgh, 'i"h"*] differs for the 
diff'erent ranges of a as characterized by by different set of [£, A, 7]. 

An acoustic black hole horizon is formed at the inner acoustic horizon as well as at the outer acoustic horizon, 
and an acoustic white hole is formed at the shock location, ~ Kin/^out can be calculated by obtaining the value 
of Kin at r)j" and Kout at i"h"*- However, it is observed that the value of Kout is of the order of 10^ times less than that 
of Kin- Also the variation of Kout is quite insensitive on a. This happens because the outer acoustic horizon forms at 
a very large distance from r+ and the Kerr parameter is essentially an inner boundary condition. Hence the 'kio — a' 
profile is almost identical with the scaled version of the 'Kin — a' profile. 

In figure 6, we plot Kin as a function of a for three different ranges of a for flow characterized by same [£, 7] 
but for three different values of A = 2.6 (uppermost panel), A = 2.17 (mid panel), and A = 2.01 (lowermost panel). 
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FIG. 5: The location of the the outer acoustic horizon rj^"* (upper panel), the stationary shock rgh (mid panel), and inner 
acoustic horizon rjf (lowermost panel) as a function of black hole spin for multi-transonic shocked flow characterized by 
[£ = 1.00001, A = 2.6, 7 = 1.43]. 

respectively. As observed for the mono transonic accretion, the surface gravity exhibits a maximum. For the multi- 
transonic flow, however, the Km^'^' at the inner acoustic horizon is characterized by a very large value of the black 
hole spin. 

In flgure 7, we plot Kin as a function of a for retrograde flow characterized by [£ = 1.00001, A = 3.3,7 = 1-4] ■ Kin 
non linearly anti correlates with a for the retrograde flow. 



B. On Possible Spectral Signature of Relativistic Acoustic Geometry 



In our attempt to make a close connection between the salient features of the acoustic geometry and any astrophys- 
ically relevant observables, we study the relation between the analogue surface gravity with shock related accretion 
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FIG. 6: The surface gravity Kin evaluated at the inner acoustic horizon rj," has been plotted as a function of black hole spin for 
three different ranges of the Kerr parameters. Multi-transonic shocked flow for a fixed set of values of [£ — 1.00001, 7 — 1.43] 
has been explored for three different values of the angular momentum A = 2.6 (the uppermost panel), A = 2.17 (mid panel) 
and A — 2.01 (lowermost panel) respectively, to cover the different ranges of the Kerr parameter. 

variables since such variables play a crucial role in determining the spectral signature of the astrophysical black hole 
candidates. Shock formation phenomena is a discontinuous event somewhat equivalent to the first order phase tran- 
sition. Accretion variables will change discontinuously at the shock location. Such accretion variables like density, 
velocity and the flow temperature provides the characteristic profiles of the observed spectra [s^l ■ A shocked flow will 
provide distinctively different spectra in comparison to its shock free counterpart. Corresponding spectra for shocked 
flow exhibits additional rich features due to the presence of the shock. The ratio of the pre and the post shock flow 
variables will determine such complex features. The fiux distribution for such spectra for axially symmetric accretion 
in the Kerr metric can be calculated in terms of such ratios (Das & Huang, in preparation). In subsequent sections, we 
will study the dependence of the values of the acoustic surface gravity on pre(post) to post(pre) shock ratio of various 
relevant quantities responsible to characterize the spectral feature - namely, the shock strength Rm = which 
is the pre to post shock Mach number ratio, the ratio of the post to pre shock density Rp = p+ /p_ (the compression 
ratio) and the temperature Rt = T+/T^, respectively, for both the prograde and the retrograde flow. In flgure 8, we 
plot the ratio of the acoustic surface gravity kio = Kin/Kout (scaled by 10~^) with the shock strength Rm = M^/M+ 
for three different ranges of the Kerr parameter for the prograde flow (top left and right and the bottom left panel) 
and a particular range of the Kerr parameter for the retrograde flow (bottom right panel). The range of the Kerr 
parameters and the values of [£,A,7] used are the same as has been used to produce the figure 6 and the figure 7. 
For low to moderately high values of a for the prograde fiow, Kjo anti-correlates with Rm in general. However, for the 
substantially large value of ^io ^ 

' profile is different. For such range of a, Kjo initially correlates with Rm and 
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FIG. 7: The surface gravity Kin evaluated at the inner acoustic horizon for multi-transonic shocked retrograde accretion 
characterized by [£ = 1.00001, A = 3.3, 7 = 1.4]. 

attains a maximum, then starts falling for larger Rm non linearly as a is increased further. For the retrograde flow, 
Kio usually anti-correlates with Rm- 

Multi-transonic flow with shocks formed closer to the black hole are associated with higher values of the shock 
strength and larger compression ratio, and hence the higher value of Rt as well. Hence the 'kjo — Rm^ profile should 
be similar with the '/tio — Rp^ and 'kio — i?T' profile. That this intuitively obvious conclusions follows in reality may 
further be formally demonstrated through figure 9 and figure 10 where we plot Kio as a function of Rp (figure 9) and 
Rt (figure 10) as well for the same set of values of the initial boundary conditions as has been used to produce figure 
8. The overall conclusion is that except for certain span of a corresponding to the near extremally rotating black 
holes, the ratio of the acoustic surface gravity at the inner and the outer acoustic horizon has relatively large value 
for weaker shocks formed in the prograde accretion. For retrograde accretion, however, the ratio Kio assumes higher 
value for weak shocks in general. 

Since the black hole spin anti-correlates with Vsh and correlates with r™, shocks formed at the larger distance 
leads to the formation of the inner acoustic horizon relatively closer to the black hole where the value of the acoustic 
surface gravity becomes substantially large. On the other hand, shocks formed closer to the black hole are rather 
strong shocks and thus produce greater compression ratio as well as the larger value of Rt- Hence weaker shocks 
correspond to the inner acoustic horizons located relatively closer to the black hole where the acoustic surface gravity 
assumes a higher value. On the other hand, as explained earlier, change of the location of the outer acoustic horizon is 
less sensitive to the infiuence of the black hole spin, hence 'Kjo — [RM,Rp, Rt]' variation is effectively equivalent with 
the scaled down version of 'k^ — [Rm, Rp, RtV- This explain why Kio anti-correlates with Rm, Rp and Rt in general. 

However, it is not quite clear at this stage why such 'k^o [Rm - Rp, RtY profile exhibits opposite nature for some 

range of the black hole spin for prograde accretion onto nearly extremely rotating holes. The inherent complexity of 
the dependence of the value of k on various quantities as evaluated on the horizon prohibits to make any conclusive 
remark in this case since the analytical expression for the quantities involved are not available within the framework 
of the formalism presented here. 

C. Dependence of k on [£',A,7] 

Both the inner and the outer acoustic horizon anti-correlates with [5, A, 7]. The variation of r™ is quite insensitive 
to £, moderately sensitive to 7 and highly sensitive to A. On the other hand, the variation of r^"* is hardly sensitive 
to A, whereas it is moderately sensitive to 7 and very much sensitive to £. The shock location rgh correlates with A 
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FIG. 8: Kio, the ratio of the surface gravities as evaluated at the inner acoustic horizon to that as evaluated at the outer 
acoustic horizon has been plotted (scaled by a factor of 10~^) as a function of the shock strength Rm for three different ranges 
of the black hole spin for the prograde flow as characterized by a fixed set of [£ — 1.00001, 7 = 1.43] and three different values 
of A = 2.6 (top left panel), A = 2.17 (top right panel), A = 2.01 (bottom left panel), and for retrograde flow as characterized 
by [£ = 1.00001, A = 3.3,7 = 1-4] (bottom right panel). 

but anti-correlates with [£,7]. As usual, stronger shocks forms closer to the black holes. 

In figure 11, the surface gravity evaluated at the inner acoustic horizon for multi-transonic shocked accretion 
is plotted as a function of £ (uppermost panel for flow characterized by [A = 2.17,7 = 1.43, a = 0.881049812]), A 
(mid panel, for flow characterized by [£ ~ 1.000004,7= 1.43, a = 0.881049812]) and 7 (lowermost panel, for flow 
characterized by [£ = 1.000004, A = 2.28, a = 0.881049812]). Kj„ anti-correlates with [£, A, 7]. The variation of Km is 
insensitive to £, moderately sensitive to 7 and considerably sensitive to A. Ko«t also anti-correlates with [£, A, 7] in 
general (not shown in the figure). Kout, however, is around 10^ to 10^ smaller in magnitude to Kin- The variation of 
Kout, however, is insensitive to A, moderately sensitive to 7 and quite sensitive to £. In the context to the acoustic 
geometry, A thus somewhat controls the properties of the flow as well as that of the acoustic geometry close to the 
actual black hole, whereas £ influences the related quantities mainly at the larger distance. 

VII. DISCUSSION 

In this work we make attempt to understand the influence of the general relativistic background black hole spacetime 
metric on the embedded perturbative manifold, i.e., on the relativistic acoustic geometry. To accomplish that task, we 
consider a specific example of classical analogue model - the equatorial slice of the general relativistic hydrodynamic 
inviscid axisymmetric accretion onto a spinning astrophysical black hole. We found that such accretion configuration 
gives rise to a very interesting example of acoustic geometry where more than one acoustic horizons may form. We 
then calculate the acoustic surface gravity for such horizon(s) and studied how the value of such surface gravity, which 
is the characteristic feature of the perturbative manifold, depends on the spin angular momentum of the astrophysical 
black holes, the Kerr parameter, which is the main characteristic feature of the background space-time - as well as 
on various accretion parameters which are the characteristic features of the background fluid continuum. Visser and 
Weinfurtner [TSj have demonstrated that the equatorial slice of the Kerr geometry is equivalent to certain analogue 
model based on a vortex geometry. Our work, in some sense, is complementary to such finding. We, however, for the 
first time in the literature provided a semi analytical formalism where the infiuence of background black hole space- 
time has explicitly been demonstrated on the embedded relativistic acous tic geom etry. Since almost all astrophysical 
black holes are supposed to posses some degree of intrinsic rotation [2^ l45l - l47l l54i [66l [69l. ItoI. IttI . the effect of the Kerr 
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FIG. 9: Kio, the ratio of the surface gravities as evaluated at the inner acoustic horizon to that as evaluated at the outer acoustic 
horizon has been plotted (scaled by a factor of 10~^) as a function of the shock strength Rp for three different ranges of the 
black hole spin for the prograde flow as characterized by a fixed set of [calE = 1.00001, 7 = 1.43] and three different values of 
A = 2.6 (top left panel), A = 2.17 (top right panel), A = 2.01 (bottom left panel), and for retrograde flow as characterized by 
[£ = 1.00001, A = 3.3, 7 = 1.4] (bottom right panel). 



parameter on the classical analogue model is extremely important to understand. Our work precisely accomplished 
that task. 

One of the most interesting aspects of acoustic surface gravity is the existence of the associated analogue Hawking 
radiation. Acoustic horizon emits Hawking type radiation of thermal phonons. The corresponding analogue Hawking 
temperature Tah = K/27r is the characteristic temperature of such radiation measured by an observer at infinity. 
Acoustic horizons explored in our work, and its associated analogue surface gravity are thus characterized by certain 
Tah- In this work, we have demonstrated how the measure of the acoustic surface gravity may be associated with 
certain observables through the spectral feature of the astrophysical black holes. This work, thus, makes the first ever 
attempt to obtain the observational signature of the analogue hawking temperature profile for a large scale relativistic 
classical fluid. One, however, should note that Tah for an astrophysical black hole being too low, might be masked 
by the background thermal noise, hence the accreting primordial black holes might serve as better candidates for the 
purpose of the possible measurement of the analogue Hawking temperature. 

The surface gravity can not be computed at the shock location since u and Cg changes discontinuously at the 
shock. As a result, (du/dr)^ ^ and (dcg/dr)^ ^ diverges. The surface gravity and the associated analogue Hawking 
temperature thus becomes formally infinite at the shock. Had it been the case that viscosity and other dissipative 
effects would be included in the system, such discontinuities in (du/dr)^ ^ and {dcg/dr)^ ^ would have been smeared 
out. Under that circumstances, k, as well as Tah would have a finite but extremely large value at the shock. This 
may be considered as a realistic manifestation of the general results obtained by jS^j . 

In our work, however, we consider only the inviscid flow. The effect of the viscous transport of the angular 
momentum, however, is quite a subtle issue in considering the analogue effects in black hole accretion. Thirty nine 
years after the discovery of the standard accretion disc theory (56i . i67i] , realistic modelling of viscous transonic accretion 
fiow with dissipation is still quite an arduous task. From the analogue point of view, viscosity is likely to destroy 
the Lorenz invariance, and the assumptions behind constructing an acoustic geometry may not be quite consistent 
for such case. Nevertheless, very large radial velocity close to the black hole implies that the infall time scale is 
substantially small compared to the viscous time scale. Lar ge r adial velocity even at larger distances are due to the 
fact that the rotational energy of the fluid is relatively low [gI l32l [soj . Our assumption of inviscid flow is not unjustifed 
from the astrophysical point of view. One of the significant effects of inclusion of the viscosity would be the reduction 
of the angular momentum. As we demonstrate in this work, the location of the acoustic horizon anti-correlates with 
A. Weakly rotating fiow makes the dynamical velocity gradient steeper leading to the conclusion that for viscous 
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FIG. 10: Kio, the ratio of the surface gravities as evaluated at the inner acoustic horizon to that as evaluated at the outer 
acoustic horizon has been plotted (scaled by a factor of 10~^) as a function of the shock strength Rt for three different ranges of 
the black hole spin for the prograde flow as characterized by a fixed set of [calE = 1.00001, 7 = 1.43] and three different values 
of A = 2.6 (top left panel), A = 2.17 (top right panel), A = 2.01 (bottom left panel), and for retrograde flow as characterized 
by [£ = 1.00001, A = 3.3,7 = 1-4] (bottom right panel). 

flow the acoustic horizons will be pushed further out from the black hole and the flow would become supersonic 
at a larger distance for the same set of other initial boundary conditions. The value of the surface gravity (and 
the associated analogue Hawking temperature) anti-correlates with the location of the acoustic horizon. A viscous 
transonic accretion disc is thus expected to produce lower value of k and Tah compared to its inviscid counterpart. 

An axially symmetric accretion flow in vertical equilibrium has been studied in our work where the disc height is 
a function of the radial distance on the equatorial plane. Axisymmetric accretion may also be studied for two other 
different flow configurations, namely, for disc with constant thickness and flow in conical equilibrium where the local 
flow thickness to the radial distance on the equatorial plane remains constant. For these two flow configurations, the 
critical points are isomorphic to the sonic points jH, [52|. Hence the location of the acoustic horizon and the values of 
the dynamical velocity along with the sound speed and their space derivative as well can be computed analytically 
without taking recourse to the integral flow solution, hence by avoiding the process of numerical integration. The 
non monotonic behaviour of the 'k — a' proflle, especially the appearance of Kmax'''^' and the dependence of of the 
corresponding amax'''^' on [£, A, 7] will then be better understood since we can then directly differentiate the expression 
of K with respect to a a to find out for what value of a (for a fixed set of [£, A, 7] ) the K — a profile attains its maximum, 
and how such value of a changes with the variation of [£, A, 7]. We plan to explore such analytical dependence in our 
future work for axisymmetric flow in conical equilibrium. 

Finally, we would like to emphasize that in the present work we did not aim to provide a formalism using which the 
phonon field generated in the system concerned could be quantized. To accomplish that task, one has to demonstrate 
that the effective action for the acoustic perturbation is equivalent to a field theoretical action in curved space, and 
the associated commutation and the dispersion relations should directly follow [T^I & [? ]. Such considerations are 
rather involved and are beyond the scope of this paper. Our main motivation was rather to apply the analogy to 
describe the classical perturbation of the fluid flow in terms of a field satisfying the wave equation in an effective 
geometry and to study the relevant consequences. 
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FIG. 11: For multi-transonic shocked accretion flow, the surface gravity Kin evaluated at the inner acoustic horizon as a 
function of the specific energy of the flow £ (uppermost panel) for flow characterized by [A — 2.17, 7 = 1.43, a = 0.881049812] , 
the specific angular momentum A (mid panel) for fiow characterized by [£ = 1.000004, 7 — 1.43, a — 0.881049812], and of the 
adiabatic index of the fiow 7 (lowermost panel) for flow characterized by [£ = 1.000004, A = 2.28, a — 0.881049812] . 



Acknowledgments 

HP and IM would like to acknowledge the kind hospitality provided by HRI, Allahabad, India, under a visiting 
students research programme. TKD would like to acknowledge the professional support provided by S. N. Bose 
National Centre for Basic Science (by offering a long term sabbatical visiting professor position) where part of the 
work has been done. The research of HP and HC is partially supported by the Natioanl Science Council of the 
Republic of China under the grant NSC 99-2112-M-007-017-MY3. The research of TKD is partially supported by the 
astrophysics project under the XI th plan at HRI. 



[1] Abraham, H., Bilic, N., Das, T. K., 2006, Classical and Quantum Gravity, 23, 2371 

[2] Abramowicz, M. A., Zurek, W. H. 1981, ApJ, 246, 314 

[3] Afshordi, N., Paczyhski, B., 2003, ApJ, 592, 354 

[4] Barai, P., Das, T. K., Wiita, P. J., 2004, ApJ, 613, L49 

[5] Barcelo, C., Liberati, S., and Visser, M., 2005, 'Analogue Gravity', Living Reviews in Relativity, Vol. 8, no. 12, websource 

http://relativity.livingreviews.org/Articles/lrr-2005-12/ , also at gr-qc/0505065| 
[6] Beloborodov, A. M., & lUarionov, A. F. 1991, MNRAS, 323, 167 






21 



[7] Bercelo, C, Liberati, S., Sonego, S., & Visser, M., 2004, New Journal of Physics, Volume 6, Issue 1, pp. 186 
[8] Bhattacharjee, J. K., Ray, A. K., 2007, ApJ, 668, 409 

[9] Bhattacharjee, J. K., Bhattacharya, A., Das, T. K., & Ray, A. K., 2009, to appear in MNRAS, see also arXiv: 0812 .4793* ^1 
[astro-ph] . 

[10] Bilic, N. 1999, Class. Quant. Grav. 16, 3953 

[11] Borozdin, K., Revnivtsev, M., Trudolyubov, S., Shrader, C, & Titarchuk, L., 1999, ApJ, 517, 367 
[12] Boyer, R. H., & Lindquist, R. W. 1967, J. Math. Phys. 8, 265 



[13] Cardoso, V., 2005, Acoustic black holes, |physics/0503042 
[14] Chakrabarti, S. K., 1989, ApJ, 347, 365 
[15] Chakrabarti, S. K., Titarchuk, L. C, 1995, ApJ, 455, 623 
[16] Chaudhury, S., Ray, A. K., Das, T. K., 2006, MNRAS, 373, 146 
[17] Das, T. K., & Chakrabarti, S. K. 1999, Classical and Quantum Gravity, 16, 3879 
[18] Das, T. K., 2002, ApJ, 577, 880 
[19] Das, T. K. 2004, Class. Quant. Grav. 21, 5253 
[20] Das, T. K., Pendharkar, J. K., Mitra, S., 2003, ApJ, 592, 1078 
[21] Das, T. K., Bilic, N., & Dasgupta, S. 2007, JCAP, 6, 9 

[22] Das, T. K., & Czerny, B., 2012, New Astronomy, Volume 17, Issue 3, p. 254-271 

[23] Dasgupta, S., Bilic, N., and Das, T. K., 2005, General Relativity & Gravitation (GRG), Volume 37, Number 11, 1877-1890 
[24] Dauser, T., Wilms, J., Reynolds, C. S., & Brenneman, L. W., 2010, MNRAS, 409, 1534 
[25] Gammie, C, & Popham, R., 1998, ApJ, 498, 313 

[26] Frank, J., King, A., Raine, D., 2002, Accretion Power in Astrophysics, Cambridge University Press, Cambridge 
[27] Fukue, J., 1983, PAS J, 35, 355 
[28] Fukue, J., 1987, PAS J, 39, 309 

[29] Goswami, S., Khan, S. N., Ray, A. K., Das, T. K., 2007, MNRAS, 378, 1407 
[30] Hawking, S. W., 1974, Nature, Volume 248, Issue 5443, pp. 30-31 

[31] Hawking, S. W., 1975, Communications In Mathematical Physics, Volume 43, Issue 3, pp. 199-220 
[32] Igumenshchev, I. V., & Beloborodov, A. M., 1997, MNRAS, 284, 767 
[33] Kafatos, M., Yang, R. X., 1994, MNRAS, 268, 925 

[34] Kato, S., Fukue, J., & Mineshige, S., 1998, Black Hole Accretion Disc, Kyoto University Press. 
[35] Kinoshita, S., Sendouda, Y., & Takahashi, K., 2004, Phys. Rev. D., 70, 123006 
[36] Laurent, P., & Titarchuk, L., 1999, ApJ, 511, 289 
[37] Liang, E. P. T., Thomson, K. A., 1980, ApJ, 240, 271 

[38] Liberati, S., Sonego, S., & Visser, M., 2000, Classical and Quantum Gravity, Volume 17, Issue 15, pp. 2903-2923 

[39] Lu, J. F. 1985, A & A, 148, 176 

[40] Lu, J. F. 1986, Gen. Rel. Grav. 18, 45L 

[41] Lu, J. F., Yu, K. N., Yuan, F., Young, E. C. M., 1997, A & A, 321, 665 
[42] Lu, J. F., & Gu, W. M., 2004, Chin. Phys. Lett., 21, 2551 

[43] Mach, P., 2009, Reports on Mathematical Physics, vol. 64, issue 1-2, pp. 257-269 
[44] Mach, P., & Malec, E., 2008, Physical Review D, vol. 78, Issue 12, id. 124016 
[45] Martinez-Sansigre, A., & Rawlings, S., 2011, MNRAS, 414, 1937 

[46] McClintock, J. E., Narayan, R., Davis, S. W., Gou, L., Kulkarni, A., Orosz, J. A., Penna, R. F., Remillard, R. A., & 

Steiner, J. F., 2011, Classical and Quantum Gravity, Volume 28, Issue 11, pp. 114009 
[47] Miller, J. M., Reynolds, C. S., Fabian, A. C, Miniutti, G., & Gallo, L. C, 2009, ApJ, 697, 900 
[48] Moncrief, V. 1980, ApJ. 235, 1038 

[49] Muchotrzeb, B., & Paczynski, B., 1982, Acta Actron. 32, 1 

[50] Muchotrzeb, B., 1983, Acta Astron. 33, 79 

[51] Muchotrzeb-Czerny, B., 1986, Acta Astronomica, 36, 1 

[52] Nag, S., Acharya, S., Ray, A. K., Das, T. K., 2012, New Astronomy, Volume 17, Issue 3, p. 285-295 

[53] Naskar, T., Chakravarty, N., Bhattacharjee, J. K., Ray, A. K., 2007, Physical Review D, vol. 76, Issue 12, id. 123002 

[54] Nixon, C. J., Cossins, P. J., King, A. R., & Pringle, J. E., 2011, MNRAS 412, 1591 

[55] Novello, Visser & Volovik (ed.) 2002, Artificial Black Holes. World Scientific, Singapore. 

[56] Novikov, I., & Thorne, K. S. 1973, in Black Holes,eds. c. De Witt and B. De Witt(Gordon and Breach, New York). 

[57] Pariev, V. I., 1996, MNRAS, 283, 1264 

[58] Peitz, J., Appl, S., 1997, MNRAS, 286, 681 

[59] Proga, D., & Begelman, M. C, 2003, ApJ, 582, 69 

[60] Ray, A. K., 2003a, MNRAS, 344, 83 

[61] Ray, A. K., 2003b, MNRAS, 344, 1085 

[62] Ray, A. K., Bhattacharjee, J. K., 2002, Phys. Rev. E, 66, 066303 

[63] Ray, A. K., Bhattacharjee, J. K., 2005a, A dynamical systems approach to a thm accretion disc and its time- dependent 

behaviour on large length scales, 2005, arXiv: astro-ph/OSllOlSvl 
[64] Ray, A. K., Bhattacharjee, J. K., 2005b, The Astrophysical Journal, Volume 627, Issue 1, pp. 368-375. 
[65] Ray, A. K., Bhattacharjee, J. K., 2006, Indian Journal of Physics, 80, 1123, also at eprint arXiv : astro-ph/ 0703301 1 
[66] Reynolds, C. S, Brenneman, L. W., Lohfink, A. M., Trippe, M. L., Miller, J. M., Reis, R. C, Nowak, M. A., & Fabian, 

A. C, 2011, 'Probing Relativistic Astrophysics Around SMBHs: The Suzaku AGN Spin Survey', arXiv: 1112 . 0036v l 



22 

[astro-ph.HE] . 

[67] Shakura, N. I., Sunyaev, R. A., 1973, A&A, 24, 337 
[68] Shrader, C, & Titarchuk, L., 1998, 499, LSI 

[69] Tchekhovskoy, A., Narayan, R., & McKinney, J. C, 2010, ApJ, 711, 50 

[70] Tchekhovskoy, A., & McKinney, J. C, 2012, 'Prograde and retrograde Black Holes: Whose Jet is More Powerful?', 

arXiv:1201.4385vl [astro-ph.HE] 
[71] Unruh, W. G. 1981, Phys. Rev. Lett. 46, 1351 
[72] Unruh, W. G. 1995, Phys. Rev. D. 51, 2827 

[73] Unruh, W. G. & Schiitzhold, R. 2005, Physical Review D, vol. 71, Issue 2, id. 024028 

[74] Visser, M. 1998, Class. Quant. Grav. 15, 1767 

[75] Visser, M. & Weinfurtner, S., 2005, Class. Quant. Grav. 22, 2493 

[76] Yang, R. X., Kafatos, M., 1995, A&A, 295, 238 

[77] Ziolkowski, J., 2010, Memorie della Societ Astronomica Italiana, v. 81, p. 294 

[78] Hereafter the phrases 'acoustic' and 'analogue' will be used synonymously for the sake of brevity. 

[79] It is interesting to note that from the historical viewpoint, however, the first ever study on the analogue system was perhaps 

accomplished by Moncrief (1980) [3 which was essentially performed on the accreting astrophysical fiows. 
[80] A positive Bernoulli's constant fiow is essential to study the accretion phenomena so that it can incorporate the accretion 

driven outflows ( [l3| and references therein). 
[81] The effect of gravity is determined within the full general relativistic framework only up to several gravitational radii. 

Beyond a certain length-scale it asymptotically follows the Newtonian regime. 
[82] For a full classification of the integral curves passing through various categories of sonic points, see, e.g.. Das & Czerny 

2012. 



